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RESUMEN 

I En base a datos observacionales para una muestra de galaxias dominadas por 

\ materia obscura y un par de ciimulos de galaxias, hemos encontrado que la densidad 

, central del halo no depende de su masa y el radio del niicleo es proporcional a la 

velocidad de rotacion maxima. Estas leyes de escala estan de acuerdo con halos 
CDM cuyas partes internas se expandieron por inestabilidades gravotermicas si las 
partfculas autointeractiian eficazmente solo en estas partes. Encontramos que la 
seccion de choque es inversamente proporcional a la dispersion de velocidades. 

ABSTRACT 

. ^ Using observational data for a sample of dark matter dominated galaxies and 

■ two cluster of galaxies, we have found that the central halo density does not depend 

\ on its mass, and the core radius is roughly proportional to the maximum rotation 

velocity. A good agreement with these scaling laws is obtained for CDM halos 
I whose dense inner parts were expanded by gravothermal instabilities if the particles 

i efficiently self- interact only in these parts. We find that the particle cross-section 

\ is inversely proportional to the velocity dispersion. 
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1. INTRODUCTION 



During several decades, dynamical studies of galaxies, and group and cluster of galaxies have pointed out 
to the existence of massive dark matter halos. On the other hand, according to current models of structure 
' formation in the universe, luminous galaxies should form from the gas trapped within the deep gravitational 
wells of dark matter (DM) structures emerged from primordial density fluctuations. In these models, non- 
dissipative, cold, colhsionless particles (cold dark matter, CDM) were required. The CDM structure formation 
scenario succesfully accounted for a wide range of observations, in particular on large scales. However, on small 
scales, compared with observations this scenario seems to predict too centrally concentrated halos and too 
much substructure in Milky Way-size halos. These discrepancies have induced to introduce some modifications 
to the CDM scenario, in particular, to the nature of the DM. 

From the point of view of particle physics, a large list of candidate DM particles has been proposed but 
unfortunately, none of the particles that might constitute the universe's missing mass have been detected 
nowadays. Nevertheless, it is possible that astronomical observations may help us to constrict some of the 
properties of these particles. For example, as was mentioned above, the existence of soft cores in the dark 
halos appears to be not compatible with colhsionless CDM particles. Therefore, astronomical studies about the 
halo properties — in particular of their cores — are crucial for understanding the nature of the dark particles 
and the structure formation in the universe, as was emphasized in a pioneering paper on this subject by J. 
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Kormcndy (1990; sec also Kormendy 1988). Here, we summarize the halo core scaling relationships we have 
inferred from observations from dwarf galaxies through cluster of galaxies (Firmani et al. 2000a, b), and we 
dicuss some of the implications of our results on the nature of the DM particles and the formation of halos. 

2. HALO CORE SCALING LAWS FROM OBSERVATIONS 

Analysis of the original virialized halo mass distribution for most of galaxies is uncertain due to the ambigu- 
ities in the estimate of the stellar mass-to-light ratios M/ L and the gravitational pull the collapsing gas exerted 
over the inner parts of the halo. This is why we limited our study to only galaxies (i) strongly dominated by 
DM and (ii) with accurately measured rotation curves. The sample taken from the literature consists of six 
dwarf galaxies, nine LSB galaxies, and two late- type low luminosity galaxies (Firmani et al. 2000b). In all 
these cases, the galaxies are DM dominated. Even so, we have substracted from the observed rotation curve 
the small disk contribution. The halo components were fitted to a non-singular isothemal model; thus, for 
each galaxy characterized by its maximum circular velocity Vmax we estimate its central density pc and core 
radius rc- For the less DM dominated LSB and low-luminosity galaxies, we have roughly calculated the factor 
by which the halo component was "deformed" due to the disk pull over the DM using the adiabatic invariance 
approximation (see details in Firmani et al. 2000b). On galaxy cluster scale, we have used the surface mass 
distribution for the cluster CL0024-I-1654 derived from unprecedent high-resolution strong lensing mass maps 
(Tyson, Kochanski, & DeU'Antonio 1998), and for the cluster CL0016-I-16 derived from weak lensing studies 
(Smail et al. 1995). In both cases there is no evidence of a massive cD galaxy and the inner mass distribution 
is soft. 

In Figs. 1 we show the dependence of on V„iax we have found from the observational data. Although 
with a big scatter, within a large range in Vmax we estimate that: 

Pc{r) w 0.02 Mgpc"^ and Tc » 5.5 

Similar results were found for an uniform sample of high and LSB galaxies of the Coma Ursa Mayor cluster 
(Verheijen 1997; §6). In this case, the rotation curve decompositions were made assuming M/Lk constant 
for all galaxies, and the halo component was fitted to a pseudo-isothermal model. In contrast, from a sample 
of Sc-Im and dwarf galaxies, Kormendy (1988,1990) inferred that pc decreases with the galaxy luminosity (or 
Vmax)- Certainly, more efforts should be done in the future in order to increase the sample of objects and to 
reduce the uncertainties in the rotation curve decomposition techniques. We remark the importance of strong 
gravitational lensing studies in order to directly probe the inner regions of the cluster of galaxies. 

3. IMPLICATIONS OF THE INFERRED HALO CORE SCALING LAWS 

The existence of soft halo cores and even more, the scaling laws obtained for DM dominated systems 

[eq.(l)], are in complete disagreement with the predictions of CDM models. Warm dark matter (WDM) has 
been proposed in order to solve the other conflict of the CDM scenario — the overlying number of guest (satellite) 
halos in a Milky Way-size halo. Cosmological N-body simulations have shown that the latter problem is indeed 
solved for a filtering scale in the power spectrum of ^ 0.1 Mpc which corresponds to a warm particle of ^ 1 
KeV (Colfn, Avila- Reese, & Valenzuela 2000). These authors have also shown that the density profiles of halos 
with masses much larger than that corresponding to the filtering scale (~ 10^ h~^M0) are very similar to those 
of the CDM models (sec also Moore ct al. 1999). Thus, even if the halos with masses near or smaller than the 
filtering mass would have a core, the more massive WDM halos will not obey the scaling laws inferred from 
observations in §2 (see also Avila- Reese, Firmani, & Hernandez 1998). 

Spergel & Steinhardt (2000) suggested other modification to the nature of the DM particles: the introduction 
of self-interaction. In Firmani et al. 2000a, the gravothermal expansion was proposed as the mechanism able to 
produce soft cores in self-interacting CDM halos. The inner velocity dispersion profile of these halos raises with 
radius. Therefore, if particles are collisional, heat transfers inwards, the core expands and cools, exacerbating 
even more the temperature gradient. This process is similar to the postcollapse gravothermal oscillations in 
globular clusters (Bettwieser & Sugimoto 1984; Goodman 1987). For globular clusters, the core expansion 
halts when the inner dispersion velocity profile flattens; this occurs because there is also an outwards heat flux 
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Fig. 1. The core radius Rc vs. Vmax,ioo for dwarfs (filled squares), LSB galaxies (empty squares) and the 
CL0024+1654 cluster (filled circle). 



from the maximum of the velocity profile. In the case of DM halos, we propose a collisional cross-section a 
such that self-interaction is efhcient only in the more dense halo regions. Besides, as the soft core grows, the 
core density decreases and at some moment, self-interaction should become inefhcient even in the inner regions. 
On the other hand, it is important to beer in mind that the CDM halo does not form instantaneously, but 
by a hierarchical mass aggregation process which establishes a cuspy inner structure with a positive velocity 
dispersion gradient. 

Recent numerical simulations for a halo with Hernquist density profile and with relatively small cross- 
sections per unit of the particle mass mx (c* = cr/nix) have shown that the gravothermal processes act in 
time scales that depend on the value of a^, (Burkert 2000; Kochanek & White 2000; see also Quinlan 1996). 
An important constriction is that the halo lifetime should be in between the core expansion time and the core 
collapse time; otherwise either the shallow core still have not been formed or the core is already in its collapse 
phase. In Firmani et al. (2000a, b), using the average observed pc and supossing that the collision time tcoi 
in the core is close to the Hubble time, a lower limit for cr^, was estimated^: cr* « 4 10^^^V|,7ax loo cm^/GeV, 
where Vmax,ioo is Vmax in units of 100 km/s. An important point to be noted is that ct* depends on Vmax or 
the maximum velocity dispersion, i.e. the cross-section is a function of the particle energy as in other classical 
physical interactions. For velocity dispersions corresponding to galaxy clusters, this value is close to the limit 
estimated by Miralda-Escude (2000) from the observationally inferred ellipticity of the cluster MS2137-23. 

One may think that the evolution of the coUisionless DM halo occurs in scales of dynamical times, tdyn, 
while those central regions of the halo affected by the gravothermal processes, evolve in relaxation time scales, 
irci- The final halo density profile is the result of both dynamical processes. The simulations carried out 
by Burkert (2000) and Kochanek & White (2000; KWOO) are for a halo already virialized. Therefore, these 
simulations do not describe the cosmological process of halo collapse and virialization. Kochanek & White 
find that the gravothermal core collapse occurs in scale times less than ~ 5 times the core formation time tc 
independent of the value of ct*. On the other hand, tc oc l/cr*. Thus, for cr* small enough the halos may still 
be in their core expansion phase. Besides, if <7* depends on Vmax as we have inferred from observations, then 
larger halos should be today in earlier stages of gravothermal expansion than smaller halos, i.e. their central 
densities have not decreased too much. This, combined with the fact that in the hierarchical scenario smaller 
halos are intrinsically more concentrated than larger ones, could produce the invariance of pc with the halo 
scale. 

^Here we assume that Vmax ~ frms.max; in fact, for CDM halos Vmax is roughly 1.3-1.7 times larger than frms.max 
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The simulation of KWOO are for a Hernquist halo, and they express cr* in unities of rfj/Mh, where Mh is the 

halo mass and rn is the scale radius of the Hernquist profile. Fitting the Hernquist profile to halos obtained 

in a N-body CDM simulation, one finds that r'j^/Mh is roughly constant. In order to obtain more quantitative 

estimates, we have used results for a ACDM0.3, h=0.7 model (Avila-Reese et al. 1999). We calculate rn as 

fH = Tv/ ch, where the virial radius r„ is defined as the radius where the average halo density is Ac times the 

background density (for our cosmology, Ac = 340), and ch is the concentration parameter which ultimately 

depends on the halo mass or Vmax and is the only free parameter in the cosmological halo density profiles. 

From the results of the simulation, we find on the average Ch = 37.5/(Vi„ax/kms~^)°'^^. The virial radius is 

1/2 

proportional to V^,, the circular velocity at this radius, and for the Hernquist profile, V^, — 2Vi„axC/^ /{1 + Ch)- 
We find that r]j/Mh « 7 lir"^^ cmVGeV. Thus, in KWOO cr* would be fi- x 7 IQ-^-* cmVGeV. For <7 = 1, 
KWOO find that tc ~ 1.7 dynamical times; after this the halo suffers the gravothcrmal core collapse. For a 
value of (T* as that we have inferred from observations, for a Vmax ~ 100 km/s halo for example (cr* « 4 10~^^ 
cm^/GeV), a ~ 0.05. The dynamical time (as deffined in KWOO) for a Vmax ~ 100 km/s halo is ^ 5 10^ years. 
Therefore, the core formation time would be of the order of a Hubble time. Halos larger than Vmax ~ 100 
km/s would have even larger core formation times. 



4. CONCLUSIONS 



• The halo core scaling laws inferred from observations of dwarf galaxies to galaxy clusters show that 
does not depend on the halo mass or Vmax and the core radius is roughly proportional to Vmax- 

• If the dark particles are self-interacting with not very large cross sections, then gravothcrmal processes 
may produce a soft core in the DM halos. Using the observational data, we estimated the value of cr* and found 
that is roughly proportional to Vj^^ax c< ^^mL.max- 

• Results from numerical simulations of already virialized halos with self-interaction, show that if ct^, is of 
the order we inferred from observations, then tc for small halos is close to the Hubble time, while for larger 
halos, tc is probably even larger, i.e. these halos are still in early stages of gravothcrmal expansion. Numerical 
simulations and theoretical studies of collapsing and virializing DM halos where self-interacion is efficient only 
in the more dense inner regions are necessary in order to attain more quantitative conclusions. 
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